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We report on a search for the flavor-changing neutral current decays D° — > e + e~ and D° — » fi + fj,~ , 
and the lepton-flavor violating decay D° — > e ± /X^. The measurement is based on 122 fb -1 of data 
collected by the BABAR detector at the PEP-II asymmetric e + e~ collider. No evidence is found for 
any of the decays. The upper limits on the branching fractions, at the 90 % confidence level, are 
1.2 x 10~ 6 for D° e+e", 1.3 x 10~ 6 for D° -> /i + fi~ , and 8.1 x 10" 7 for D° -> e ± ^ . 

PACS numbers: 13.20.Fc,11.30.Hv,12.15.Mm,12.60.-i 



In the Standard Model (SM), the flavor-changing 
neutral current (FCNC) decays D° — > e + e~ and 
D° — > 0] are highly suppressed by the Glashow- 

Iliopoulos-Maiani (GIM) mechanism Q- Their decay 
branching fractions have been estimated to be less than 
10~ 13 even with long-distance processes included. This 
prediction is orders of magnitude beyond the reach of 
current experiments. Furthermore, the lepton-flavor vio- 
lating (LFV) decay D° — » e ± fi T is strictly forbidden in 
the SM 0. 

Some extensions to the Standard Model can enhance 
the FCNC processes by many orders of magnitude. For 
example, i?-parity violating supersymmetry can increase 
the branching fractions of D° — > e + e~ and D° — > 
to as high as 10~ 10 and 10 -6 , respectively [4j. The same 
model also predicts the D° — > e ± /i T branching fraction 
to be of the order of 10~ 6 . The upper bounds on the 
predicted branching fractions of D° —> and D a —> 

e ffl are close to the current experimental sensitivities. 
As a result, searching for the FCNC and LFV decays in 
the charm sector is a potential way to test the SM and 
explore new physics. Similar arguments hold for rare K 
and B decays, but the charm decay is unique since it is 
sensitive to new physics coupling to the up-quark sector. 

In this paper, we present a search for the decays of 
D° -> e+e - , D° -> and D° -> e ± ^. The analy- 

sis is based on 122 fb -1 of data collected on or near the 
T(45) resonance by the BABAR detector at the PEP-II 
asymmetric e + e~ collider. 

The BABAR detector, which is fully described in 
provides charged-particle tracking through a combination 
of a five-layer double-sided silicon micro-strip detector 
(SVT) and a 40-layer central drift chamber (DCH), both 
operating in a 1.5 T magnetic field in order to provide 
momentum measurements. The identification of charged 
kaons and pions is achieved through measurements of 
particle energy- loss (dE/dx) in the tracking system and 
Cherenkov cone angle (0 C ) in a detector of internally re- 
flected Cherenkov light. Electrons are identified primar- 
ily in a segmented CsI(Tl) electromagnetic calorimeter, 
while muons are identified by their penetration through 
the iron plates of the magnet flux return. 

The charmed mesons considered for this analysis origi- 
nate from the fragmentation of charm quarks in the con- 
tinuum e + e~ — > cc process. There is no advantage in 
including D° decays from the B mesons because of their 
higher combinatoric background. The D° — > l + £~ (£ = 



e, //) branching ratio is determined by 

B(D° ^e + r)=S(N ohs -N hg ), (1) 

where iV bs is the number of D° — > £ + £~ candidates ob- 
served, 7Vb g is the expected background and S is the sen- 
sitivity factor, defined as: 

S = B(D° ^tt+tt-)-!-^ . (2) 

Here B(D° -> tt+ti--) = (1.43 ± 0.07) x 10" 3 is the 
D° — * 7r + 7r~ branching fraction |fj, N n7r is the number 
of reconstructed D —> ir + ir~ decays, e« and e WT are the 
efficiency for the corresponding decay mode. We choose 
D — » 7r + 7r _ as the normalization mode because it is 
kinematically similar to D° — > £ + £~ and therefore many 
common systematic uncertainties cancel in the calcula- 
tion of the efficiency ratio e^ir/tee- The key to the anal- 
ysis is to reduce backgrounds as much as possible while 
maintaining a high signal efficiency. 

We first outline the general event selection require- 
ments common to all the data samples used in the anal- 
ysis and later describe tighter optimized criteria specific 
to each decay modes. A pair of oppositely charged tracks 
is selected to form a D° candidate. They are fit to a com- 
mon vertex and only the candidates with fit probability 
larger than 1 % are retained. Since charmed mesons from 
e + e~ — » cc events are produced with momenta higher on 
average than those from e + e~ — > bb events, a minimum 
value of 2.4GeV/c is imposed on the center-of-mass mo- 
mentum of each D° candidate. In order to further reduce 
the background, the D° candidate is required to be from 
a D* + — > D°ir + decay. The D° candidate and the pion 
from the D* + are fit to a common vertex with a beam 
spot constraint. The probability for this fit is required to 
exceed 1 %. The resolution of the mass difference between 
the reconstructed D* + and D° candidates is approxi- 
mately 0.25 MeV/c 2 . We require that \Sm\ < 2.0 MeV/c 2 , 
where 6m = m{D°^+) - m(D°) - 145.4 MeV/c 2 . In ad- 
dition, all the tracks are required to have a minimum 
number of measurement points in the SVT and the DCH. 

We require the electron and muon candidates to have 
momenta larger than 0.5GeV/c and l.OGeV/c in the 
laboratory frame respectively. In this range, the aver- 
age electron and muon efficiencies are about 95 % and 
60 %, and their hadron misidentification probabilities are 
measured from t decay control samples to be around 
0.2% and 2.0%. Pion identification is also applied to 
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the daughters of D° — > tt + it~ decays. The correspond- 
ing single pion identification efficiency is around 90 %. No 
particle identification (P1D) is applied on the soft pion 
from the D* + decay. 

Except for particle identification, the selection criteria 
applied to the D° — > ir + ir~ mode are the same as those 
used for the D° — > £ + £~ modes. The signal efficiencies 
of D° — » £ + £~ and D° — > ir + ir~ are evaluated using 
a Monte Carlo (MC) simulation. We use PYTHIA Q 
for the fragmentation of the produced cc. The final 
state radiative effects are simulated for all decays using 
PHOTOS js. The detector response is simulated with 
GEANT4 |9j, and the simulated events are then recon- 
structed in the same manner as the data. 

Due to large final-state radiation and bremsstrahlung 
backgrounds, the invariant mass distributions, mu, of 
D° — > e + e~ and D° — * e ± ^ =F have a low mass tail. 
We define an asymmetric signal mass window (1.8045 < 
mu < 1.8845 GeV/c 2 ) for all three decay modes. The 
lower boundary of the signal mass window is chosen to 
include the majority of the radiative tail of D° — > e + e~ 
and D° — > e ± /i =F . The higher boundary corresponds to 
a little more than 2a of the D° mass resolution mea- 
sured from the D° — > 7r + 7r~ control sample. In order to 
avoid any possibility of bias, a blind analysis technique 
has been adopted. All events inside the D° mass window 
were hidden from inspection until the final event selection 
criteria were established and all systematic uncertainties 
were determined. 

The D° — > £ + £~ background can be taken as a sum 
of two components: a peaking background from D° — > 
hrh~ (h = K,tt) decays and a combinatoric background 
from other sources. The copious two body hadronic 
D° — > h + h~ decays will mimic the dilepton signals if 
both hadrons are misidentified as leptons. MC studies 
show that only the decay D° — > 7r + 7r~ contributes in the 
signal window. The D° -> #-71-+ and D° -> 
backgrounds peak in the lower mass region because of 
the high kaon mass. To estimate the number of peaking 
background events, N{£, we apply the selection criteria 
for D° -> £+£- to MC simulated D° -> tt+tt" events 
with lepton misidentification rates measured from a con- 
trol sample. 

MC studies show that the combinatoric background in 
both the signal mass window and high mass sideband 
region (1.9045 < m« < 2.0545 GeV/c 2 ) is dominated 
by the combination of two random leptons. The invari- 
ant mass distribution of the random lepton pair is flat. 
This is indeed consistent with what is observed in the 
high mass sideband of the data. As a result, the ex- 
pected combinatoric background in the signal window is 
just the number of dilepton events in the high mass side- 
band scaled by the ratio of the width of the signal region 
to the high mass sideband region. 

To further reduce the background, we added a selec- 
tion on the proper decay time, ct, of the D° candidate, 



and tightened our selections on the signal mass window 
and 8m. We determine the optimal selection criteria by 
maximizing the value e«/-/V sons , where N scns is the av- 
eraged 90 % confidence level upper limit on the number 
of observed signal events that would be obtained by an 
ensemble of experiments with the expected background 
and no real signal [ljj. Studies show that the correla- 
tions among the optimized discriminating variables are 
negligible. 

The expected combinatoric background there- 
fore can be factorized as 

lilacs RSmRct, (3) 

where N$b is the number of high mass sideband events 
passing the loose event selection criteria; i? maS s is the 
expected background rejection factor for a given signal 
mass window; R$ m and R c t are the expected background 
rejection factors for the tighter 5m requirement and the 
ct requirement respectively. In order to avoid possible 
bias due to the statistical fluctuation in the high mass 
sideband, we determine the discriminating variable dis- 
tribution shape of the combinatoric background from the 
MC and D — » tt + it~ control sample. The information is 
subsequently used to predict the background changes as 
a function of a particular set of selection criteria rather 
than by directly examining the data in the high mass 
sideband. The optimized final selection criteria are sum- 
marized in Table [i] The estimated numbers of back- 
ground events are listed in Table [H] The proper time 



Mode m ee [GeV/c 2 ] 5m [MeV/c 2 ] ct 
ee 1.8245 <m ee < 1.8845 \5m\ < 0.6 ct > 
fj.fi 1.8445 <m MM < 1.8845 \5m\ < 0.6 
en 1.8445 <m^< 1.8845 \5m\ < 0.5 



TABLE I: The summary of the optimized event selection cri- 
teria of D° -> i+r. 

requirement is found to be useful only for the ee mode. 
The background in the ee mode is dominated by combi- 
natorials with zero average lifetime and is halved, with a 
reduction of less than 20 % in signal efficiency, by requir- 
ing that the proper time of the D° candidate be positive. 
Such a requirement is not applied to the fifi mode, where 
its reflection background has large contribution and a 
similar lifetime behavior to the D° meson, or to the e/i 
mode, where the background is very small. 

As an important check of the background estimate, we 
have compared the expected distribution in the low mass 
sideband ( 1.6545 < m u < 1.8045 GeV/c 2 ) with the 
data. The peaking background in the low mass sideband 
is evaluated from the D° — ► h + h~ MC sample using pre- 
cise measurement of the lepton misidentification proba- 
bilities. The random lepton pairs are inferred from the 
events in the high mass sideband. Unlike the upper side, 
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the combinatoric background in the low mass sideband 
has contributions from the combination of two hadrons 
and the combination of one real lepton with one hadron, 
where the hadrons are misidentified as leptons. We esti- 
mate those backgrounds using MC data and known lep- 
ton misidcntification rates. We find that the predicted 
background distributions and levels in the low mass side- 
band (before and after the optimization of our event se- 
lection criteria) have excellent agreement with our obser- 
vation in the data for all three decay modes. 

The number of D° — > tt + it~ candidates in the data, 
Nnn, is extracted by fitting their invariant mass distri- 
bution with a binned maximum likelihood fit. The signal 
distribution is modeled as a double Gaussian, and the 
background distribution is approximated as a linear func- 
tion. The number of reconstructed D° mesons is found 
to be between 7000 and 12000, depending on the selec- 
tion criteria. The relative uncertainties in N„ n are about 
1%. 

The invariant mass distribution of the dilepton candi- 
dates after applying the optimized event selection crite- 
ria is shown in Fig. ^ The number of events observed 
(-^Vobs) and the expected background (iVb g ) are shown in 
Table [H] with no significant excess found in any decay 
mode. 

The largest systematic uncertainty in the signal effi- 
ciency ratio e^/eu calculation is due to the PID effi- 
ciency. It ranges from 1.2% for the ee mode to 4.2% 
for the fi/i mode relative to their efficiency ratio. Other 
sources of systematic uncertainty are found to be small, 
including the track reconstruction efficiencies, track mo- 
mentum resolution and MC statistics. 

The systematic uncertainties of the background esti- 
mate arise predominantly from the finite data available 
in the high mass sideband for the ee and e/z modes. For 
the im mode, a large fraction of the background is pro- 
duced by misidentified D° — > ir + n~ decays. The rel- 
ative uncertainty associated with the estimate of muon 
misidcntification is found to be about 4.7%. 
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0.02 
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0.21 
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iV bg 


2.23 ± 0.38 


4.63 ± 0.45 


2.14 ±0.36 


S [1(T 7 ] 


2.25 ±0.12 


4.53 ±0.30 


3.27 ±0.20 


N obs 
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1 





UL obtained 


1.2 x 10 -6 


1.3 x 10 -6 


8.1 x 10 -7 



TABLE II: The summary of the number of expected back- 
ground events (iVb g ), the sensitivity factor (S), number of ob- 
served events (Nobs), and the branching fraction upper limits 
at the 90 % confidence level for each decay modes. The uncer- 
tainties quoted here are total uncertainties. The uncertainty 
of iV^g is negligible for the ee and e/i decay modes. 
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FIG. 1: The dilepton invariant mass distribution for each 
decay mode. The dashed lines indicate the optimized signal 
mass window. 



The branching fraction upper limits (UL) have been 
calculated including all uncertainties using an extended 
version of the Feldman-Cousins method [ljj. All of 
the uncertainties have a negligible effect on the limits. 
The results are listed in Table ILTI 

In summary, we have performed a search for the FCNC 
decays D° — ► e+e - , D° — ► and the LFV decays 

D° — ► e+zi^ using the BABAR detector. No evidence is 
found for these decays. The upper limits on the branch- 
ing fractions at the 90 % confidence level are 1.2 x 10 -6 for 
D° -> e+e', 1.3 x 10~ 6 for D° -> and 8.1 x 10~ 7 

for D° — > e /x^. Our result improves the present best 
limits by a factor of 5 for the ee mode a little less 
than 2 for the [ifi mode ^jj, and 10 for the e/z mode [lij . 
The upper limits for the branching fractions of the e/i 
and yit/x modes begin to confine the allowed parameter 
space of R-parity violating supersymmetric models . 
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